Vortex-lattice melting in a one-dimensional optical lattice 
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We investigate quantum fluctuations of a vortex lattice in a one- dimensional optical lattice. Our 
method gives full access to all the modes of the vortex lattice and we discuss in particular the Bloch 
bands of the Tkachenko modes. Because of the small number of particles in the pancake Bose- 
Einstein condensates at every site of the optical lattice, finite-size effects become very important. 
Therefore, the fluctuations in the vortex positions are inhomogeneous and the melting of the lattice 
occurs from the outside inwards. Tunneling between neighbouring pancakes substantially reduces 
the inhomogeneity as well as the size of the fluctuations. 

PACS numbers: 03.75.Lm, 32.80.Pj, 67.40.-w, 67.40.Vs 



Introduction. — Very rapidly rotating ultracold 
bosonic gases have been predicted to form highly- 
correlated quantum states l2Lbj|- I n these states, the 
Bose-Einstein condensate has been completely depleted 
by quantum fluctuations, and quantum liquids appear 
with excitations that can carry fractional statistics. Some 
of these states have been identified with (bosonic) frac- 
tional quantum Hall states 0, 0, • It is by now a long 
standing goal to observe the experimental signatures of 
these very interesting states in the context of ultracold 
quantum gases. The conditions for these states to form 
have been expressed in the requirement that the ratio 
v = N/N v of the number of atoms N and the number of 
vortices N v , should be smaller than some critical value 
v c . The ratio v plays the role of the filling factor and 
estimates for the critical v c are typically around 8 0, Q • 
These estimates are made for infinitely large systems. 

However, observed filling factors are up till now always 
greater than 100, where almost perfect hexagonal lattices 
form and no sign of melting can be seen (jj. These ex- 
periments are carried out with Bose-Einstein condensates 
consisting of typically 10 5 particles, whereas the maxi- 
mum number of vortices observed is around 300. De- 
creasing the number of particles results in loss of experi- 
mental signal, whereas the number of vortices is limited 
by the rotation frequency that has to be smaller than 
the transverse trapping frequency. Adding a quartic po- 
tential, which stabilizes the condensate also for rotation 
frequencies higher than the transverse trap frequency, 
has until now not improved this situation although 
it has opened up the possibility of forming a giant vor- 
tex in the center of the cloud In recent work, how- 
ever, it was realized that quantum fluctuations of the 
vortices can be greatly enhanced, without loosing exper- 
imental signal, by using a one-dimensional optical lattice 
[HI- The optical lattice divides the Bose-Einstein con- 
densate into a stack of two-dimensional pancake conden- 
sates that are weakly coupled by tunneling as schemat- 
ically shown in Fig. Q The number of particles in a 
single pancake is much smaller than in a Bose-Einstein 
condensate in a harmonic trap, and therefore the fluctu- 



ations are much greater. Moreover, by varying the cou- 
pling between the pancakes, it is possible to study the di- 
mensional crossover between decoupled two-dimensional 
melting, and the strong-coupling limit where the three- 
dimensional situation is recovered. Recently, the density 
profiles for quantum Hall liquids in such a geometry have 
also been calculated [H| . Because of the small number of 
particles in each pancake shaped Bose-Einstein conden- 
sate, finite-size effects become very pronounced in this 
setup. In particular, the critical filling factor for the 
melting of the lattice v c changes compared to the ho- 
mogeneous situation. Moreover, melting is not expected 
to occur homogeneously but starts at the outside and 
then gradually moves inwards as the rotation speed in- 
creases. Therefore, phase coexistence is expected, where 
a vortex crystal is surrounded by a vortex liquid. In this 
Letter we study this interesting physics by investigating 
the quantum fluctuations of the vortex lattice for realis- 
tic numbers of particles and vortices. To decide whether 
or not the vortex lattice is melted we use the Lindemann 
criterion, which in this case has to be applied locally. 

Besides the possibility of observing highly-correlated 
quantum states after the vortex lattice has melted, also 
the vortices themselves have very interesting properties. 
For moderately rotating gases, the vortices are observed 
to order in a Abrikosov lattice 0, E| . The lowest mode 





FIG. 1: (Color online) Setup in which the melting of the vor- 
tex lattice is studied. An optical lattice along the z-direction 
(indicated in blue) divides the condensate into pancake con- 
densate which are coupled by tunneling processes. 
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of the lattice is the so-called Tkachenko mode U£\. which 
recently has been investigated experimentally 0, [l5[ and 
theoretically 0, 0, ^| . Since the study of the fluctu- 
ations of the vortex lattice requires full information on 
all the modes of the lattice, we also can identify this 
mode. Because of the optical lattice, this mode actually 
becomes a Bloch band and we also determine the disper- 
sion of the Tkachenko modes in the axial direction of the 
optical lattice. 

Lattice vibrations. — The action describing the system 
in the rotating frame is given by S = J dt J <i 3 x £(x, £), 
with the Lagrange density given by 



^!-y(x) + ^-||^| 2 



where ^(x, t) is the Bose-Einstein condensate wavefunc- 
tion, m is the mass of the atoms, Q is the rotation fre- 
quency, L z = ih(yd x — xdy) is the angular momentum 
operator, and g = 4ith 2 a/m is the interaction strength, 
with a the three-dimensional scattering length. The com- 
bination of an optical lattice in the axial direction and 
harmonic confinement in the radial direction gives an ex- 
ternal potential V(x) = V z cos 2 (Az/2) + \muj\{x 2 + y 2 ), 
where A is the wavelength of the laser producing the op- 
tical lattice. Assuming that the optical lattice is deep 
enough we can perform a tight-binding approximation 
and write ^(x,t) = $^ #(z - Zi)S&i(x,y,t), where 
is chosen to be the lowest Wannier function of the optical 
lattice. In this way we obtain the Lagrange density 
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where (ij) denotes that the sum is taken over nearest 
neighbouring sites. We have defined the interaction co- 
efficient g' = 4:7rh 2 a/\ /r 27r£ z m and hopping amplitude 
t = 4Vz /4 El /4 expl-2-y/Vz/Ez], where E z = 2ir 2 h 2 /\ 2 m 
is the recoil energy associated with the optical lattice and 
l z = X/2tt(V z /E z ) 1 ^ 4z is the harmonic length associated 
with the optical lattice. 

Melting is only expected for a Bose-Einstein conden- 
sate that is weakly interacting. The wavefunction can 
then be taken to be part of the lowest Landau level. 
That is to say we consider wavefunctions of the form 
Y\i(z - Zi)exp[-\z\ 2 /2], where z = (x + iy)/£, Zi = 
(xi + iy%)/£ and (xi,yi) is the position of the i th vor- 
tex. Here £ is the 'magnetic length'. To increase the 
validity of this study, we use i as a variational param- 
eter instead of fixed it to the radial harmonic length, 
such that our results are also valid for stronger inter- 
actions The associated frequency is uj = hjmi 2 . 
From now on distances are rescaled by i, frequencies 
are scaled by the radial trappig frequency cu±, and we 
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FIG. 2: (Color online) Classical vortex lattice and density 
profile for rotation frequencies Q/oo± = .8, .9, .95 and .97. 
Here U = 10, which corresponds to N = 250. White means 
high density, black low density. The vortex positions are in- 
dicated by a red dot, such that also the vortices outside the 
condensate are visible. 



define a dimensionless interaction strength by means of 
U = Ng' /47r£ 2 hj± = Na/V2ir£ z . In the calculations 
we take the scattering length of 87 Rb, A = 700 nm and 
V z /E z = 16, which gives U = 25N. The on-site en- 
ergy part of the Lagrange density depends in this ap- 
proximation only on the particle density and becomes 
(u + l/u- 2n)r 2 n(r)/2 + 2ttuUti 2 (y). The lowest Lan- 
dau level wavefunctions and, therefore, also the atomic 
density, are fully determined by the location of the vor- 
tices. To consider the quantum mechanics of the vortex 
lattice we, therefore, replace the functional integral over 
the condensate wavefunctions by a path integral over the 
vortex positions. 

Next we want to determine the quadratic fluctuations 
around the Abricosov lattice. To do so, we first have to 
find the classical groundstate. We calculated this ground- 
state for up to 37 vortices. For small numbers of vortices, 
the groundstate is distorted from the hexagonal lattice 
[2l| . In general, there are also vortices far outside the 
condensate. The coarse-grained atomic density is well 
approximated by a Thomas- Fermi profile [2^. For fixed 
interaction U and different rotation frequencies pictures 
of the classical groundstate are given in Fig [21 Then we 
study the quadratic fluctuations by expanding the action 
up to second order in the fluctuations 0. This yields 
an action of the form 



S = ^Ui.(Tm-E) 
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where = (5xi ti , 6y 2 ,i, Sx 2i u *2/2,i, • • •) = (. . . ,u n », . . .) 
is the total displacement vector of all the point vortices 
on site z, and T, E and J are matrices depending on f2, 
Z7, and the classical lattice positions. To diagonalize this 
action along the z-axis, we perform a Fourier transfor- 
mation to obtain 

S = u k ' ( Tidt - E - t(l - cos[fcA/2]J) • u k (4) 

fe 

Finally, we completely diagonalize this action by the 
transformation v& = P^u^ such that the action becomes 
S = J2k,a v ka( id t-v a (k))v ka where u a (k) are the mode 
frequencies of the vortex lattice. This means that the 
Vkof) where k labels the momentum in the ^-direction 
and a labels the mode, correspond to bosonic operators 
with commutation relation [^ a ,% a ^] = fc'ta'- This 
allows us to calculate the expectation value for the fluctu- 
ations in the vortex positions, but also for the correlations 
between the various point vortices. 

Tkachenko modes. — The Tkachenko modes are 
purely transverse modes of the vortex lattice. In a har- 
monic trap with spherical symmetry they become modes 
which are purely angular. In the radial direction their 
spectrum is discretized, because of the finite lattice size. 
The number of radial Tkachenko modes equals the num- 
ber of vortex rings. For 37 vortices 6 Tkachenko modes 
can be identified. A close comparison with continuum 
theory for a finite-size system, where also a discrete spec- 
trum was found [2^, is possible but beyond the scope 
of this Letter. Moreover, the Tkachenko modes also 
have a dispersion in the ^-direction. Without the op- 
tical lattice some aspects of these modes were recently 
investigated [24[. For typical parameters this disper- 
sion is plotted in Fig. [3 As is clearly visible, there 
is one gapless mode, which is linear at long wavelengths, 
while the other modes are roughly just tight-binding like. 
Moreover, various avoided crossings between these modes 
are clearly visible. The gapless mode is the Goldstone 
mode associated with the spontaneously broken rota- 
tional 0(2)-symmetry due to the presence of the vor- 
tex lattice. When the tunneling rate is very small, the 
gapped modes have exactly a tight-binding dispersion 
and the gapless mode gets a dispersion proportional to 
sin [k A/2]. This can be understood by observing that in 
this case the modes are decoupled and the Hamiltonian 
for the gapless mode reduces to the Josephson Hamilto- 
nian H = -E c d 2 /d$ + £j $>os(& - ^) 2 , which 
indeed has this dispersion upon quadratic expansion. 

It is interesting to note that for a small rotation 
frequency, which implies a small vortex lattice, the 
Tkachenko modes are not the lowest-lying modes. For 
U = 10, a Tkachenko mode becomes the lowest-lying 
gapped mode when Q > .978, but there are many modes 
in between the second and the third Tkachenko mode. 
This confirms the expectation that increasing the vortex 
lattice will bring down the Tkachenko spectrum more and 
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FIG. 3: Dispersion of the Tkachenko modes along the direc- 
tion of the optical lattice. For this plot the parameters where 
chosen as: U = 10, t = 1/10, and Q = .97. One gapless 
linear mode and five gapped tight-binding like modes can be 
identified. The gapless mode is the acoustic Goldstone mode 
associated with the broken 0(2)-symmetry due to the pres- 
ence of the vortex lattice. 



more. 

Vortex- lattice melting. — Quantum fluctuations of the 
vortices ultimately result in melting of the vortex lattice. 
To decide wether or not the lattice is melted, we use 
the Lindemann criterion, which in this inhomogeneous 
situation has to be applied locally. The Lindemann cri- 
terion means that the lattice is melted, when the ratio 
(u^i) / of the displacement fluctuations of vortex n at 
site i and the square of the average distance to the neig- 
bouring vortices A n i exceeds a critical value c 2 L , which is 
known as the Lindemann parameter. For easy compar- 
ison with previous work on this topic, we use the value 
c\ = .02. Because the coarse-grained particle density 
decreases with the distance to the origin, vortices on 
the outside are already melted, while the inner part of 
the crystal remains solid. Therefore a crystal phase in 
the inside coexists with a liquid phase on the outside. 
In Fig. 0] we compute the radius of the crystal phase 
R cr normalized to the condensate radius i?, as a function 
of the rotation frequency for fixed numbers of particles 
and an interaction strength U ', and for various hopping 
strengths t. Here we define the condensate radius R as 
the radius for with the angularly averaged density drops 
below .003. The crystal radius R cr is defined as the ra- 
dius of the innermost vortex ring that is melted accord- 
ing to the Lindemann criterion. When according to this 
definition R cr > R, we set the crystal radius equal to 
the condensate radius, i.e., R cr = R. The ratio R cr /R 
shows discrete steps because of the ring-like structure in 
which the vortices order themselves. We compare this 
with a simple local density calculation, where the crite- 
rion N/N v = n(r)/n v (r) = 8 derived in Ref. |(| is applied 
locally, using a Thomas-Fermi density profile to describe 
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FIG. 4: Crystal radius R cr normalized to the condensate ra- 
dius R as a function of the rotation frequency for N = 300, 
U = 10. The full and dashed line are for t = and 1/100, 
respectively. The dotted line is the result of a local-density 
approximation. In the inset the number of vortices N v that is 
within the condensate is plotted as a function of the rotation 
frequency. 

the coarse-grained atomic density. From the on-site en- 
ergy we extract that the variational parameter £ is set 
by the condition uj = O, and the Thomas- Fermi radius is 
given by R 4 = 16fl 2 U /(1 — Vt 2 ). When taking a constant 
vortex density n v (r) = 1/tt£ 2 , the crystal radius normal- 
ized to the Thomas Fermi radius becomes R C t/Rty — 

^l-m 2 /N = y/l - 16/N^n 2 U/(l - n 2 ^. For com- 
parison this line is plotted in Fig. EI As expected for 
a finite system, the melting occurs considerably earlier 
than predicted by the local-density theory. 

When the tunneling between pancakes is turned on, 
the fluctuations will also be coupled in the axial direc- 
tion. This decreases the fluctuations in the vortex dis- 
placements because the stiffness of the vortices increases, 
and melting occurs for higher rotation frequencies, as is 
visible in Fig. 0] To determine the presence of crystalline 
order in the axial direction, we calculate (^e lCL '^ UNi ~ Un ^)^ 
which is related to the structure factor. For the cen- 
tral vortex we always find long-range order, while for 
the other vortices, we obtain an algebraic decay. This is 
in agreement with the expectation for a one-dimensional 
system at zero temperature. Note that the quantitative 
results in Fig. 21 depend on the value of the Lindemann 
parameter. Changing this value shifts the curves, but the 
qualitative features remain the same. 

The liquid surrounding the crystal has rotational sym- 
metry. Experimentally this can be observed, by noting 
that in the liquid the vortices are no longer individually 
visible, while in the crystal they are still visible, although 
their positions are broadened by the fluctuations. The- 
oretically it is a challenging problem to describe the co- 
existing crystal-liquid. This will allow to decide on the 



occurrence of melting based on energy considerations and 
thus shed more light on the accuracy of the application 
of the Lindemann criterion in this inhomogeneous situa- 
tion. In the future we want to make a closer connection 
with continuum theory. Finally we want to extend the 
analysis to include a quartic potential, in which case giant 
vortex formation is prediction. In this case the liquid can 
also form in the center, instead of at the outer edge. Fi- 
nally we also want to investigated the melting at nonzero 
temperatures, which is experimentally relevant, because 
the zero-temperature limit is difficult to reach |25j |. 
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